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Abstract 

We derive Bezout identities for the minimal polynomials of a finite sequence and use 
them to prove a theorem of Wang and Massey on binary sequences with a perfect 
linear complexity profile. We give a new proof of Rueppel's conjecture and simplify 
Dai's original proof. We obtain short proofs of results of Niederreiter relating the 
linear complexity of a sequence s and K(s), which was defined using continued 
fractions. We give an upper bound for the sum of the linear complexities of any 
sequence. This bound is tight for sequences with a perfect linear complexity profile 
and we apply it to characterise these sequences in two new ways. 

Keywords Bezout identity, continued fractions, linear complexity, minimal polynomial. 

1 Introduction 

In [11] we showed how to obtain a minimal polynomial (MP) of a finite sequence recur- 
sively. We begin with a compact matrix reformulation of this: multiplying an updating 
matrix and the current 'MP matrix' gives the new one. The product rule for determinants 
gives a quick proof of Bezout identities for MP's (which was proved from first principles 
in pT|). Next we give several basic characterisations of sequences with a perfect linear 
complexity profile (PLCP) in terms of MP's. 

Section 14.21 applies the Bezout identities to give a new proof of a theorem of Wang 
and Massey characterising binary sequences with a PLCP. 

We give a new proof of Rueppel's conjecture and simplify the proof of p[]. This 
result is that the linear complexity (LC) of the first n terms of the binary sequence 
r = (1, 1, 1 , 1, 3 , 1, 7 , 1, ...) is L^nJ- We will see that the updating matrices of r are 
either a constant matrix U or an identity matrix. (When applying the Euclidean algorithm 
to r of length a power of two, the successive remainders can be obtained from the matrices 
M, U, . . . ,U, M, [U Lemma 3]; we obtain M as the MP matrix of the first two terms of 
the sequence.) Our approach is to work directly with a family of binary polynomials 
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{^,W : i > 0} rather than with the roots of Y 2 + xY + 1 in some algebraic extension 
of F 2 (x)[F], as in [I]. (In fact, the polynomial 9 n of [3j equals xy( n ' and the requisite 
lemmas of [1] can be easily deduced from ours.) This means that Rueppel's conjecture 
can be proved more simply: use the Euclidean algorithm and properties of the {7^}. 

Section 14.41 studies a quantity K(s) which was defined using the continued fraction 
expansion of the generating function of s in [7j. In [9], the author shows that K(s) equals 
the supremum of the quantities e n which figure prominently in our main theorem (giving 
MP's recursively). We take their supremum as our starting point, defining the height of s 
to be ht(s) = sup{e n }. We give short proofs of inequalities which ht(s) satisfies (originally 
proved for K(s) using continued fractions) as well a characterisation of PLCP sequences 
in terms of K(s) which appeared in [TJ. In this way, we can deduce results of [7J, [9] as 
corollaries. 

We conclude with an upper bound for the sum of the LC's of a sequence. This bound 
is tight for PLCP sequences and leads to two new properties of these sequences which are 
equivalent to Rueppel's original definition. 

We thank an anonymous referee for Theorem 14. 251 which improves an earlier approach. 
Some of the results of this paper were presented in May 2010 at Projet Secret, INRIA, 
Rocquencourt, France. We would also like to thank the project members for their interest 
and hospitality, and Nicolas Sendrier for a useful question. 

2 Preliminaries 

2.1 Notation 

The letter n always denotes n a strictly positive integer, N = {1, 2, . . .}, Z = {0, ±1, . . .}, 
F denotes a field, F q denotes a finite field of order q and D is a commutative integral 
domain with 1^0. For any set S containing 0, S x = S \ {0}. 

As usual, -D[[x _1 ,x] is the domain of Laurent series in x" 1 , D[x~ l ,x] is the domain 
of Laurent polynomials and v : — > Z U {—00} is the exponential valuation: 
v(0) = —00 and if g 7^ 0, v(g) = max{i : gi 7^ 0}, so that v coincides with deg on D[x}. 

It is elementary that v(gg') = v(gr)+v(gr'), v(g + g') < max{v(g),v(g')} and v{g + g') = 
max{v(g), v(g')} if v(g) 7^ v(g)'. We also use v denote its restriction to D\x~ x ,x\. 

2.2 Sequences 

An infinite sequence is a function s : N — > D. The set of infinite sequences over D clearly 
forms an abelian group. We can regard it as a submodule of a natural -D[x]-module 
as in P21 Section 2] as follows. We begin with as standard D[[x _1 ]-module i.e. 

acting on itself via multiplication. This also makes into a -D [a;] -module. Let 

s = J2j>i s j a; ~' J - Then for / G D[x], put 

i>i 
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One checks that o makes x 1 D[[x x ]] into a D[x]- module. 

Definition 2.1 An infinite sequence s satisfies a linear recurrence relation if it is a tor- 
sion element i.e. z/Ann(s) = {/ G D[x] : / o s = 0} ^ {0}. In other words, if for some 
f G D[x] x , (/ • s)d-j = for d — j < — 1 where d = deg(f) i.e. 

foSj-d H h /dSj = /or d + 1 < j. 

When /d = 1, we can write Sj = — (/o%-d H h fd-iSj-i) for j > d + 1 and s is a linear 

recurring sequence. 

A finite sequence is s = (s\, . . . , s n ) G -D n and 

S = S n X~ n + • • ■ + SiX^ 1 . 

For 1 < i < n, we write for (si, . . . , Sj). In the following definition, the multiplication 
of / G _D[a;] and s G a;" 1 /^^ -1 ] is in 

Definition 2.2 (Annihilator) / 7i0j Definition 2.7, Proposition 2.8]) We say that f G 
D[x] is an annihilator (or a characteristic polynomial) of s G D n if f — or {f ■ s)d-j — 
for d — n < d — j < — 1 i e. 

/oSj-d H H /dSj = /or d + 1 < j < n 

where d = deg(/) > ; written f G Ann(s). 

The following definition is a functional version of [10l Definition 2.10] A : D[x] x x 
D n+1 — > D is given by 

d 

&(f, S) = (/ • S)rf_ n _i = /fe Sn+i-^+fc 

fc=0 

where d = deg(/). If s G -D n+1 is understood, we write A n+1 (/) for A(/, s); if / is also 
understood, we simply write A„ +1 . Clearly A n+1 (l, s) = s n+ i and / G Ann(s) if and only 
if / G Ann(sW) and A n+1 = 0. 

3 MP- matrices of a sequence 

The principal result in this section is Theorem 13.51 
3.1 Minimal Polynomial 

Any polynomial of degree at least n annihilates s G D n vacuously, so Ann(s) ^ (0) and 
the following definition makes sense. 
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Definition 3.1 (Minimal Polynomial (MP)) (JJH Defn. 3.1]) We say that f G 
Ann(s) is a minimal polynomial (MP) of s G D n if 

deg(f) = min{deg(a) : g G Ann(s) x }. 

Let MP(s) denote the set of MP's of s. We do not require MP's to be monic. The linear 
complexity (LC) of s is L(s) = deg(/) where / G MP(s). We will also write L n for L(s) 
when s is understood; similarly Lj = L(s^) for 1 < j ' < n. It is convenient to set L = 0. 
The following two functions simplify many statements in what follows. 

Definition 3.2 (Exponent Function) We define the exponent function e : D[x] x x 
({0}UN) -»■ Z 

e(./»=n + l-2deg(/). 

When / is understood, we write e n for e(/, n)\ for example, we often write e n = n + 1 — 2L n 
if / G MP(s). We conventionally set eo = 1. 

Definition 3.3 (Index Function) Let s G D n . We set /i (0) = I. Given G MP(s (i) ) 
and Aj + i = A J+ i(//^) /or 1 < j < n — 1, put 0' = —1 and define n' inductively by 

, _ j n — 1 if A n 7^ and e n _x > 
L (n — 1)' otherwise. 

To simplify statements in the remainder of the paper, we will write /x'^ -1 - 1 for /z^-? -1 )') 
and A; for A^'" 1 ), ^O'- 1 )^ 1 )). 

The polynomial part of f and s, written [/ • s], is the Z)[x]-summand of / • s: 

v(s)+d 
3=0 

When s is understood, we will write [/] for [f • s\, 



3.2 MP matrices 

Definition 3.4 Let s G D n . We will call a 2 x 2 matrix M = M(s) with entries from 
D[x] an MP-matrix for s if M n G MP(s), M 12 = [M u ], M 21 = M' n and M 22 = [M 21 \. 

When s is understood, we will write for M^(s). To express the main theorem in 
matrix terms, we need an integer Heaviside function 

1 if i > 
v if i < 0. 

Trivially, e • 9(e) = max{e, 0}. 
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Theorem 3.5 (Cf. ffl, JEtf) Let s E D n , e E D be arbitrary but fixed. Put 



1 

£ -1 



and Ao = 1. Suppose that is an MP-matrix for for 1 < j < n — 1 and e = e n _i. 
If A n = then M^™ -1 ) is an MP-matrix for s, A' n+1 = A^ and e n = e + 1. 
On the other hand, if A n ^ then 
(i) L n = max{e, 0} + L n _i = n' + 1 - L n / 

zs an MP-matrix for s where 



Tj(n-l) 



A' • £ e ' e ( e ) — A n • x _e '( 1_e ^ e ^ 
™#(e) l-0(e) 



e„ = — |e| + 1 



A n if e> 
K ife<0 



v ) ifMlf ^ 0, deg(M^ J ) = max{e, 0} + deg(M- 



r(") 



(n-l)> 
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Example 3.6 Let r = (1, 1, 0, 1, 0, 0) E F% and e = 0. We have M (0) = J fry definition, 
e = 1 and A]. = ri = 1, so tnat AfW = [/^M^ = f J = [/ say. iVe:rf ; e x = and 



1 



1 1 
1 



giving 



A 2 = r 2 = I, so that U {2) = 

M (2) = U(2) M (1) 

In the same way, e 2 = 1 = A 3 and 

M (4) = M (3) = ^(2) = 

We will also write M^ 1 ' in terms of MP's as 

M (n) 



"11" 




x 1 




" x + 1 1 " 


1 




1 




1 



X + X + 1 .T 
X + 1 1 



[/i <»>] 



(n) 



>{n) [../(n 



where /z^ E MP(s) and // n ) = (/i^, [a^™- 1 ])- In this formulation, Theorem 13.51 yields 

A' n ■ x e p*"- 1 ) - A n • p 7 ^" 1 ) if e > 

A; • JI^-V - A n ■ x~ e Jjf^ otherwise 



and we see again that there is a net increase of e in LC precisely when e > 0. For the 
convenience of the reader we recall the algorithm implied by Theorem I3.5| it is a rewrite 
of [EH Algorithm 4.6]. 
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Algorithm 3.7 Algorithm MP ( Cf. % p. 184]), 13 P- 124]) 

Input: n > 1, e G D and s = (s%, . . . , s n ) G D n . 
Output: ZZ = (/!, [//]). 

{ e: =l;7l':= (e, -1); A' := 1; ZZ := (1, 0); 
FOR j = 1 TD n 

{A := Efc=o^ s fc+ 2+£; 
IF A ^ 

THEN{IF e < THEN /I := A' • /Z - A • a;- e /Z'/ 

ELSE {(*,«) := /Z; 

ZZ: = A' ■ a: e ZZ - A ■ /Z / / 
Z/:=(^,u); A':=A; 
e := -e}} 

e := e + 1} 
RETURN (ZZ)} 



The analogue of M^ n ' in Berlekamp's context and notation is 



a uj 

t 7 



[H p. 181]. See 



pO') 



also [3, p. 180]. We derive our analogue as follows. Let p : {0, . . . , n — 1} — > {0, . . . , n — 1} 
be defined by 

"1 if J = 

j — j' otherwise. 

Then p(n) = p(n— 1) + 1 if A„ = 0. Let jj^* and jj'^* denote the reciprocal polynomials. 
Considering the cases e > and e < gives the following corollary of Theorem 13.51 

Corollary 3.8 If A n ^ then 

(i) M W * = A; • i£ n -V * - A n • x^™- 1 ) n'fr-Q * 

(ii) p(n) = p(n — 1) + 1 if e„_i < and p(n) = 1 otherwise. 

In this way, we obtain a linear feedback shift-register of shortest length L = n+ \~ e and 
'feedback polynomial' generating s^ 1 '. Our corresponding updating matrix is 

A'„ -A n ■ 
9(e) l-9(e) 

where p = p(n — 1), p(0) = 0, p(n) is set to if e > and then p(n) = p + 1 (regardless 
ofe). 



3.3 Bezout Identities 

Definition 3.9 For s G D n , we set Vo = 1 and 

V n _i if A n = 



A^ +1 V„_i z/A n ^0. 
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Table 1: Algorithm MP e = 0, r = (1, 1, 0, 1, 0, 0) e F|. 



3 


A? 


e i-l 









i 




1 





1 


i 





a; 


1 


2 


i 


1 


x + 1 


1 


3 


i 





X 2 + X + 1 


X + 1 


4 





1 


X 2 + X + 1 


X + 1 


5 


i 





x 3 + x 2 + 1 


x 2 + X + 1 


6 





1 


x 3 + x 2 + 1 


x 2 + X + 1 



Proposition 3.10 (Cf. J|J Theorem 7.42]) \M^ \ = -V and \M^\ = V„. 

Proof. The first statement is trivial. Inductively, if A n = 0, there is nothing to prove; 
otherwise Theorem 13.51 gives 

| M (n)| = |jy(n-l)| . | M (»-1)| = A ' n+i . |Af(™-l)|. 



Thus we have n Bezout identities for s: for 1 < j < n 

- ^ = Vj (1) 

proved from first principles in [JTJ Theorem 3.3]. Again we see that gcd(/i^, = 1 if 

D has unique factorisation; likewise, gcd(//^, fj/V') = 1. 



4 PLCP Sequences 

The following is a slight generalisation of [T2]: s G D n has a perfect linear- complexity 
profile (PLCP) if = [^J for 1 < j < n. If s : N -> D, then s has a PLCP if for all n, 
has a PLCP. It is easy to see that the binary sequences of length 1 to 4 with a PLCP 
are (1), (l,s 2 ), (1, 1, 0), (1, 0, 1) and (1, 1, 0, s 4 ), (1, 0, 1, s 4 ); see Table SJfor their Jl^. 



4.1 Basic Characterisations 

Recall that for any sequence, fi^ = 1 and e = 

Proposition 4.1 The following are equivalent: 

(i) s has a PLCP 

(ii) Li = 1 and for 2 < j < n 



Lj — Lj_ i 



i/j is even 

1 otherwise, 



Table 2: G F 2 [x] 2 for < z < 4. 



J 







(1,0) 


1 


CM) 


2 


(x+ Ai,l) 


3 


(x 2 + Aix + l,x) 


4 


(x 2 + (A 1 + A 3 )x + l,x + A 3 ) 



(m,) /or 1 < j < n 

1 z'/j z's even 
otherwise, 



(iv) Aj j£ for all odd j, 1 < j < n 

(v) for 2 < j < n, 

u - iy 

(vi) J2^ = (x + e, 1) and /or 2 < j < n 

= 



j — 2 if j — 1 is even 
j — 3 otherwise, 



Aj_i • ^ — Aj • /i^ 2 ) z/j z's even 



Aj_ 2 ' xfi^ ^ — Aj ■ i$ 3 ** otherwise. 

Proof, (i) (ii) (iii): Easy consequence of the definitions. 

T tI = T = L i = ^'j- 1 



(i) =>• (iv): If j < n + 1 is odd then Aj 7^ 0, for otherwise ^- + l = 2^- = L j = L J 



(iv) (i): Let A, ^ for all odd 1 < j < n + 1. Then s ± ^ 0, L x = 1 and 
ei = 0. If A2 = 0, then L2 = Li = 1, otherwise L2 = max{ei, 0} + 1 = 1, so that L2 is as 
required. Suppose that j < n is odd and L& = [^"J ^ or all /c, 1 < A; < j — 1. We have 
Lj = j—Lj-i = j — 3 ~Y- = L^irJ • If J = n+1, we are done. Otherwise, if A J+1 = 0, we have 
L j+ i = Lj = |_^jrj = L^J, whereas if A J+1 ^0,L i+1 = J + 1- Lj = j + l- = L^J- 

(iii) =>- (v): Let i = (J — 1)'. We have 

If j is even then Aj_i 7^ and ej_ 2 = 1, so z = j — 2; if j is odd then ej_ 2 = and 
ej-3 = 1 so z = (j — 2)' = j — 3. (Note that z + 1 is always odd, so that A,; + i 7^ 0.) 

(v) =>■ (vi): Inserting (j — 1)' and e 3 _i in Theorem 13.51 gives the formulae for fx. 

(vi) =>- (iii): We have Li = 1, Lj = Lj-i if j is even and Lj = L 3 -_i + 1 if z is odd. Thus 
if j is odd, ej = j + 1 — 2L,- = j + 1 — 2(Lj_i + 1) = j + 1 — 2L ? _i = ej_i + 1. Applying 
this inductively gives (iii). ■ 

Thus without loss of generality, we may assume that n is odd, and it is easy to see 
that for odd n, /i^ + c • /i /(n) = /i (n) + c • ^ {n ~ l) G MP(s) for any c E D. 
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Table 3: <x (j) G F 2 [x] for < % < 4. 



3 


(70 





1 


1 


X + 1 


2 


(A! + 1)X+1 


3 


(Ai + l)x 3 + x + 1 


4 


(Ax + l)x 3 + (AiA 3 + l)x + 1 



Corollary 4.2 Let D = F be a field, s G D n , fj,^ and Aj be as in Theorem \3.5\ for 
l<j<n. If s has a PLCP and G MP(s), then 



Li) 



z/j zs even 



(x + c)- /iC?- 1 ) - . ^0-3) ^ 



i-2 



erwise 



where c E D. 

Proof. If j is even, then has a unique monic MP and if j is odd, any monic MP 
of is 



A, 



A, 



M (i) + c . M 'C?) = _ . M (i-3)) + c . M (i-D = ( x + c ) . M 0-D _ . ^0-3) 



A,-_ 2 



for some c 6 D by [TOl Theorem 4.16]. ■ 
4.2 The Characterisation of Wang and Massey 

Here we prove a theorem of Wang and Massey [H] on binary sequences using the Bezout 
identities ([I]) and results of Section |4~T1 Thus D = F 2 throughout this subsection. Let us 
call s G D" stable if s x = 1 and for even j, 2 < j < n, = sj + Si • 

The transform t which appears in the next result was used in [El Theorem 3] and a 
similar one was used in 1 141. 



Proposition 4.3 Let n be odd and t = s 2 + (x + l)s + 1. Then s is stable if and only if 
tj = for j even, < j < n. 



Proof. We have to — Si + 1 and tj = Sj + s,-+i + s L for all even j, 2 < j < n — 1. 



9 



a (j) 



Lemma 4.4 (Cf. Lemma 1]) Let a® = [p,W 2 -t). If s has a PLCP then for 2 <j <n, 

+ Aj ■ o-^'- 2 ) + Aj ■ (x + 1) if j is even 

x 2 aV~ l > + cr^~ 3 ^ + x(x + 1) otherwise. 
Proof. From the definition, 

a® = (x + + /i^ 2 + [/x&>] 2 . 



Let j be even, so that //^ = l > + A • 2 ) by Theorem 14. 1[ where A = Aj. Putting 
W = [/z^] and expanding cr^ gives 



which is as required since p^~ 2 ' jjtP" 1 ' + n^~ 2 'u^~ x ' = 1 by Proposition 13.101 The proof 
for j odd is similar. ■ 

The reader may also check that deg(ck')) = j — 1 if j is even and j if deg(cr^) = j is odd, 
but we will not need this. 

Theorem 4.5 Let n > 1 be odd and s G D n . If s has a PLCP then s is stable. 

Proof. An easy inductive proof using Lemma 14.41 shows that if s has a PLCP, then for 
< j < n, Oq = 1 and = for 2 < j < n. Let j be even with < j < n — 1. As 
to = 0, we can assume inductively that tj = for even j, 2 < j < n — 3. Put /j = fj^ n \ 
which has degree Now 

= [fi 2 ■ t] = [fi 2 ■ {t 3 x~ 3 + ■■■ + t n ^x x - n + t n x~ n )\. 

As = 0, the quadratic term in the right-hand side VIZ. jjj n + 1 ■ t n _i is zero. This forces 
t n _i = as /in+i is the leading coefficient of /i, and s is stable by Proposition 14.31 ■ 



Proposition 4.6 The number o/s6F" which have a PLCP is (q — 1) Tal qrL 2 J . 

Proof. A sequence s determines (Ai, . . . , A n ) G D n uniquely, and conversely. Thus the 
result follows from Proposition 14. l( iv) . I 



Corollary 4.7 flTQ Theorem p. 13]) Let n > 1 be odd and s G D n . Then s has PLCP if 
and only s is stable. 

Proof. As in noted in [H] . there are clearly |~|] stable sequences in D n , so the result 
follows from Proposition 14.61 ■ 

For other proofs of the characterisation of Wang and Massey, see [7J Corollary 1] 
(which uses the characterisation of [1]) and [H Theorem 3] (which uses an idea of [13]). 
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4.3 Rueppel's Conjecture 

Throughout this section, r : N — y F 2 denotes the binary sequence with r 2 k = 1 for k > 
and zero otherwise. We have already seen the two invertible matrices 



U 



x 1 
1 



and M = M (2) (r) 



x + 1 1 
1 



as well as powers of U in Example 13.61 Our goal in this section is to show that the 
updating matrix of r is U if n is even and the 2x2 identity matrix if n is odd. 

To work with the powers of U, we define a sequence of binary polynomials as follows: 

7 (o) = Qj 7 (i) = 1 and 

7 (*) = ary (*-i)+ 7 (*-2)f 0r j b >2. 

The first eight are 0, 1, x, x 2 + 1, x 3 , x 4 + x 2 + 1, x 5 + x and x 6 + x 4 + 1. A simple 
induction gives the powers of U in terms of the polynomials 



Proposition 4.8 For k > 1, 



U 



(fc+i) (fc) 
7 (fe) 7 (*-i) 



Lemma 4.9 (TJf. Lemma 1]) 

(i) For m > n, 7 ( m+n ) = x 7 ( m ) 7 (™) + 7 <m ~ ri -). In particular, 7 ( 2n ) = or/™) 2 , 
fwj deg( 7 ^^) = n — I. 

(m) ( 7 (") + 7 ("- 1 ))(o) = 1. 

(ty) gcd^"),^" 1 )) = 1. 

Proof, (i) The case n = 1 is the definition. Suppose inductively that the result is true 
for n — 1 < m. We proceed as follows: 



7 



(m+n) 



(m+n-1) _|_ (m-n-2) 



U^M/yCn-l) _|_ 7 (m-n+l)j _|_ ^(m+n-2) 

x 7 (m) [x 7 ( n -^ + 7 (n ~ 2) ] + x 7 (m) 7 ( "- 2) + x 7 (m -" +1) + 7 ( m +"- 2 ) 

X ry{m)j(n) _)_ X7 ( m ) 7 ("-2) _|_ x 7 ( m ~ n+1 ) -f 7 ( m +™- 2 ) = X7 ( m ) 7 (") _|_ 7 ( m -«) 



,(2n-2) 



7^ 



since if m = n, x^ n '^ n > + xx ' + 7 
and we are done. Otherwise the inductive hypothesis yields 



(2n-2) + (2) + (1) + 7 (2n-2) _ q _ „(m-n) 



X7 ( m ) 7 ("--2) = 7 ( m +™~2) _|_ 7 <m-n+2) _ 7 (m+n-2) _|_ (m-n+1) _|_ ^(m-n) _ 

The remaining items are easy inductions. 



Proposition 4.10 For > 0, (%) 7 ( 2 ^ = x 2 " 1 and (ii) 7 



,(2 fe -l) 



2 fc -2 J 
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Proof, (i) The result is true for k = 0. The rest of the proof is a simple induction 
using m = n = 2 k ~ l in Lemma l4.9( i). The proof of Part (ii) is similar to [5J Lemma 
1(3)]. The result is true for k = 0, so suppose inductively that it is true for k — 1. Since 
2 fe _ 2 i = 2i + • • • + 2 k -\ we have 



k-l 



k-l k-l 



k-l k-l 



J2 x 



2 fc_ 2 j 



9 fc _ 



(2*) 



We claim that the right-hand side is 7^ l \ This is true for k = 1 and assuming that it 
is true for k — 1, the right-hand side is 



1 + 



fc-l fc-2 



1 + 



(2 l ) 



(2 fc - 1 ) 



• X7 V " ' = 1 + 7 1 " ~' • X7 V 



( 2 fc-i_i) _ ( 2 fe-i) 



1 + ( 7 (2fc - 1} + 1) 



j=i i=j 

where the last equality follows from Lemma I4.9( i) 
We have 



L lo g2 n \ 

£ 

i=0 



x 



The next lemma is essentially [U Lemma 2] with a simpler proof. It is key to determining 
the discrepancies. 



Lemma 4.11 For k > 1 



,(2 fc ) n 



,t M U 



l TT 2-2 k 



1 

X + 1 



Proof. For p > 2 direct evaluation yields 



p 


1 






X + 1 





7 



(p+i) + 7 (p-i) 



Proposition 14.101 with p = 2 k implies that 



7 (2'-l) +7 (2») = I 2 



2 fe_ 2 j 2 fc 
x = x r 



and 7 (2k+1) + 7 (2fe " 1} 



(i)^(2 fe ) 



X7 V '7 



x 



Now we can state and prove the main result of this section. 



Theorem 4.12 (Cf. If MW = M^{r) then = M and for n > 3, 

M^™ -1 ) ifn is even 

C/M^- 1 ) = U^rM ifn is odd. 



M in) 
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Proof. The matrix M = M^> and the result for n = 3, 4 were derived in Example 
Suppose inductively that n is odd and the result is true for sequences of length n — 1 > 4. 
First we use Lemma [4.111 to show that /i*™ --1 ) = /i*™ --2 ) ^ Ann(r^ n ^). Let k = |_log 2 nJ so 
that r*™) = A 2 \ Lemma [4.111 and the inductive hypothesis give 



M (n_2) 



r (n) 
1 



n — 3 



lf- p x 



1 

x + 1 



a; 



7 0-i) 7 (p) 
(p) 7 ( P +i) 



where p = 2 h 



re+l 



> 0. As L n _ 2 = the discrepancy 2 ) • r( n ^)n zl i_ n is 



1 

x + 1 



(( x+ 1) 7 (p) + 7 (p-D) +7 



(^) 



which is 1 since deg(7 ( - p - ) ) = p — 1, deg( 7 ( p ^) = p — 2 and n > 1. Next we construct M^ n \ 
As L n _x = L„_ 2 = (n - l)/2, e n _i = n - 2L n _ x = 1 and M< n > = UM {n ~^ = U^M by 
Theorem 13.51 and the inductive hypothesis. 

It remains to show that fjfi^ £ Ann(r^ n+1 )). If n + 1 < 2 fc+1 then r( n+1 ) = r^ 2 ) and 
applying the first part gives 



M (n) 



1 



UM {n ~ 2) 



,(2 fc ) 1 



7 (p-2) (p-1) 

7 (p-l) j(p) 



1 

X + 1 



So the discrepancy (jjl^ ■ r} 11 ^)^.^ is ((x + l^" 1 ) + 7 (p ~ 2) ) p+1 + 7 ^+i_ n _ 1 = 0. 
If n + 1 = 2 fe+1 then 



M (n) 



r (n+l) 
1 



C/" p x" 



n-l 


1 






X + 1 





7 (P-!) ry(p) 
(p) 7 (p+l) 



x 



-n-l 



1 

X + 1 



where p= ^=±. Then (//(") • r> +1 )) n+1 _ _ = ((x + + 7 ^- 1 )) Ii±i +7 l „ 2 +1 _ n which 

2 n 2 2 n 

is zero as before. Hence p*™- 1 £ MP(r' n+1 ') and the proof is complete. ■ 



Corollary 4.13 For the sequence r, p 



(n) a 



IS 



i/n zs even 
( 7 (p) _|_ 7 (p-i) )7 (p-i)) 2 j n ^ odd 



where p = i/ n zs odd and L n = \_^-\ ■ 

Corollary 4.14 (Cf. J^) If py- n ' is as in Corollary 4-13 , its reciprocal defines a linear 
feedback shift-register of shortest length [^^J which generates 



,(n) 



Remarks 4.15 (i) Let Y 2 + xY + 1 £ F 2 (x)[F] have roots p,p~ x in some algebraic 
extension of F 2 (x) and 6 n = p n + p~ n , as in ^j. Parts 1,2 of ^ Lemma 1} easily 
imply that for n > 1, x / y ( - n \x) = 9 n (x). Lemma 4. 11 trivially implies Lemma 2] and 
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hence Lemma 3]. Further, Lemma 4] is a trivial consequence of Proposition \4~1\ 
and [4\ Lemma 5] is immediate. It follows that the main results of J^j can be proved 
using the binary polynomials 7 ( - n ) and the Euclidean algorithm, i.e. without introducing 
Y 2 + xY + 1 G F 2 (x)[F] and its roots. 

(ii) It is not hard to show that if s : N — >■ F is a linear recurring sequence and 
f G Ann(s), then f G MP(s) if and only ifgcd(f, [/]) = 1. Also, if s G F n and f G MP(s), 
then gcd(/, [/]) = 1; apply Proposition \ 3.10\ or see ITO{ Corollary 3.24]- We note here 
that r(") shows that the converse fails for finite sequences. Let 2 k < n < 2 k+1 . Then 

so that x 2k G Ann(r^ n ^) and gcd(a; 2fc , [x 2k ]) = 1, but if n is even or n + 1 < 2 k+l then 
L ra (^) = L^^J < 2 fc and x 2 is not an MP of r^ ; cf. [7| p. 230]. It would be interesting 
to know when f G Ann(s) and gcd(/, [/]) = 1 implies f G MP(s). 

(Hi) As noted in FA p. 231], the theorem of Wang and Massey also shows that r has 
a PLCP. 

4.4 The Height of a Sequence 

Let D = F be a field, s : N — > F and let {A n G F[x]} be the partial quotients in the 
continued fraction expansion of s. Then K(s) was defined in [7J p. 223] by 

K(s) = sup{deg(A n )}. 

n>l 

The next theorem was proved using inequalities satisfied by K(s) in [9], Theorem 2] 
for Frequences. (As in [7] et seq., s : N — > F is irrational if it is not a linear recurring 
sequence.) 

Theorem 4.16 If s : N — > F is irrational then 

K(s) = sup{e n }. 

As {e n } figures prominently in Theorem 13.51 we will take their maximum as our starting 
point. In general, the range of L(s) gives 1 — n < e n < n + 1. 

Definition 4.17 (Height) If s G D n , we set 

ht(s) = maxjej : 1 < j < n} 
and if s : N — > D , we put ht(s) = lim^oo ht(s^ ra ^). 

As ht(s^) < ht(s ( ^ n+1 - 1 ), the limit always exists, although it may be infinite. For 
example, if s = (0, . . . , 0) G D n , then ht(s) = n + 1, so that if s is the infinite zero 
sequence, ht(s) = oo. We have ei = 2 if s\ — and ei = 1 otherwise, so that ht(s) > 1. 

The terminology 'height' was suggested by Theorem(i) 13.51 for s G D n , L n increases 
(by e) exactly when e > 0. Thus ht(s) is the maximum of the degree jumps in s, and 
we can compute it using Algorithm MP or the Berlekamp-Massey algorithm. Since LC is 
non-decreasing, lim n _ i>00 L n exists. 

Theorem 14. 121 immediately implies the next result. 
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Corollary 4.18 (Cf. J7j) ht(r^) = 1. 



Proposition 4.19 (Cf. Q Theorem 3]) 

(i) If s e D n then ht(s^) > ej > 1 - ht(s^) for 1 < j < n. 

(ii) Ifs:N^D then ht(s) > e n > 1 - ht(s). 

Proof. The inequality e^ < ht(s^) is trivial. We prove that 1 - ht(s^) < ej by 
induction, the case n = 1 being a trivial verification. Suppose that n > 2 and the result 
is true for all sequences of length n — 1. If A n = or (A n 7^ and e n _i < 0), then 
e n = e n _i + 1 > (1 - ht^"- 1 ))) + 1 > 1 - ht(s) since ht(s( n_1 )) < ht(s). If A n ^ 0, then 
e n = -e n _i + 1 > -ht(s( n_1 )) + 1 > 1 - ht(s). Now let s : N D. If ht(s) = 00, the 
result is trivially true, otherwise taking limits gives the required result. ■ 

If s has a PLCP then ht(s) = 1 = e 2 n and Q2n-i = = 1— ht(s) so that the bounds of 
Proposition 14. 19l are tight. 

The case k — 1 in the next two results relates to [TJ Theorem 2] and p5] Corollary 1]. 

Proposition 4.20 If s E D n or s : N — > D and 

( 1 — k if n is odd 
n ~ 1 k if n is even 

for some fixed integer k > 1 then ht(s) = k. 

Proof. The definition implies that ht(s) < k and Proposition 14. 191 gives 1— ht(s) < 1 — k. 



Proposition 4.21 (Cf. /7, Theorem 2], fM, Corollary 1]) If lim^oo L n (s) = 00 and for 
all n, e n > 1 — k for some fixed integer k > 1, then ht(s) < k. 

Proof. Suppose that for some n, e n > k + 1. Let m > n be the first integer for which 
A m +i 7^ (such an m exists since lim n _ i , 00 L n = 00). Then e m > e n > k + 1, so that 
e m+i = — e m + 1 < —k which is a contradiction. Hence e n < k for all n and ht(s) < k as 
required. ■ 



Proposition 4.22 (^C/. |2l Theorem 3]) If s G £> n or s : N D is irrational then s has 
a PLCP if and only ifht(s) = 1. 

Proof. If ht(s) = 1 then by Proposition ITOl < e n < 1 i.e. f < L n < 2±1 and 

L n = L^^J • The converse follows from Theorem 14.11 ■ 



Remark 4.23 We can now use Propositions \4-19 
results of ffl, |f. 



4-22 and Theorem 4-16 to deduce 
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4.5 The LC Sum 

To simplify the notation in this section, we put Lj = Lj(s), a = and 



cr n (s 



E L - 



i=l 



where s G -D n . We will make repeated use of [61 Theorem 2]: if A n 7^ then L n = 
max{L n _x, n — L„,_i} or L n = max{e, 0} + L n _! in the notation of Theorem l3.5( i). We will 
need the fact that ^ILiL^rJ = J which is easily proved by induction. We begin 

with a technical lemma. 



Lemma 4.24 For integers k > — 1 and I > 1, 

k+2l 



E 

i=k+l 



i + 1 



l 2 + (k + l)l. 



Proof. If we put m — k + 1 + 1, the sum is 



2-1 

E 

i=0 



m — i 



+ 



m + i + 1 



2-1 

E 

i=0 



m — i m + i + 1 1 
— :— + ^ ^ 



km 



since m — i and m + i + 1 have opposite parity. 



Theorem 4.25 For s G D n , a n < Y.LA^l 

Proof. It is convenient to set <t_i = L_! = 0. Let us call j > —1 stable if it is 
odd, Lj = and Oj < J2i=il^\- Clearly —1 is stable. Suppose inductively that 

k = 2c — 1 > — lis stable, so that L^ = c. We show that (i) we can assume that there is 
an / > 2 such that Lj = c for k + 1 < i < k + / < n, Lk+i+i ^ c and either (ii) k + 21 < n 
and k + 21 is stable or (iii) n < k + 21 and n is stable. This will complete the proof as we 
can replace k by k + 21 and c by c + / in (i) until for some k and I, n < k + 21. Then (iii) 
applies and we conclude that n is stable. 

(i) By Theorem I3.5( i). L fc+1 = c independently of A fc+1 . We are done if n — k+ 1, so 
suppose that n > k + 2. If A fc+2 7^ then L k+2 = c + 1 = |_ fc+ ^ +1 j and we can replace k 
by A; + 2 and c by c + 1. So we can assume that A fc+2 = 0. Hence for some I > 2, Lj = c 
for A; < i < k + I. If k + / = n, then <r n < X]™ =1 L^J since k is stable. So we can assume 
that k + 1 < n and L k+l+1 7^ c. 

(ii) Let n > k + 21. Firstly, Afc + ; + i 7^ since L k +i+i 7^ c. Hence L^+^+i = c + / by 
Theorem I3.5( i). As c + / = max{c + /, k + 21 — (c + /)}, we have 

J c for A; < i < A; + 1 , , 

4 ~ \ c + / for fc + Z < i < k + 21 [ ' 
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and in particular, Lk+21 = j . Equation fl2]) and Lemma [4.241 imply that 



k+2l 



k+2l 



u = ic+ i(c + 1) = i 2 + (k + 1)1 = 



i=k+l 



i=k+l L 



i + 1 



and since k is stable, o~k+2i = o~k + I 2 + (k + 1)1 < Yli=i |_^"J i- e - ^ + 2/ is stable. 

(iii) Suppose now that n < k + 21 and let m = n — k — I — 1, so that < m < / — 1. 
As is stable, it is enough to show that YH=k+i < Sr=fc+i L^J • But 



fc+i— m— 1 



E L < = E L * + E( L 



k+l-i 



L 



=fe+i 



i=k+l 



i=0 



and for k+l<i<k+l— m — 1, Lj = c< |_^rJ • Further, Equation ([2]) implies that the 
second summand is (m + l)(2c + /) = (m + 1)(A; + / + 1). As in the proof of Lemma T4.24[ 



E 

i=0 



+ I - i + 1 



+ 



k+l+i+2 



(m + !)(£; + / + 1). 



Therefore ^™ =fc+1 Lj < XliLfc+i L^J an d the proof is complete. ■ 
The following consequence of Theorem 14.251 appeared in [5] for even n and D a field. 

Corollary 4.26 For s G D", J™=i U < L( n + X ) 7 4 J ■ 

Two more equivalent conditions for a PLCP follow. 

Corollary 4.27 Let s G -D n . T/ie following are equivalent: 
(1) s has a PLCP 

(ii) U < LirJ for 1 < % < n and a n = [^^\ 
(m) U > forl<i<n. 

Proof. Theorem 14.251 shows that (i) =>- (ii). (ii) =>- (i): The case n = 1 is a trivial 
verification. Suppose inductively that the converse is true for sequences of length n—1 > 1, 
U < L^J for 1 < i < n and a n = L^^J • Then [^J + L^J 



Or, 



°~n-l + L„ < 

L^J + L n by Theorem WM We conclude that L„ = [^J and so a n ^ = J . The 
inductive hypothesis now shows that s^" 1 ) has a PLCP. 

(i) =>- (iii) is trivial, (iii) =>- (i): This is trivial for n — 1. Suppose inductively that 
the converse is true for sequences of length n — 1 > 1 and that L« > L^jrJ f° r 1 <i <n. 



Then a n > IXiL^J = ^-^^l and therefore a n = [ (J ^\ by Theorem MB By the 



(n+l) 2 



same argument, cr n _i = |_ij-J anc ^ hence L n = [ 4 j — |_irJ = L 2"J - This together with 
the inductive hypothesis implies that s has a PLCP. ■ 
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Example 4.28 Let s be the geometric sequence (1, 1, 1) with minimal polynomial x + 1 
and Li = L 2 = L 3 = 1. As L 3 ^ 2, s does not have a PLCP. In fact both (i) a 3 = L~r^J 
and (ii) L3 > 2 fail. Now lett = (1, 1, 1, 0) which does not have a PLCP since s does not. 
Algorithm MP gives x 3 + x 2 + 1 G MP(t) ; so L 4 = 3. We have a 4 = 6 = L^r^-J , but both 
(i) L4 < 2 and (ii) L 3 > 2 fail. 



Remark 4.29 Let D be a field F. Corollary easily gives an upper bound for the 
number of ¥ -multiplications required by Algorithm MP on a sequence of length n. We can 
clearly divide by A^ +1 at iteration i, which requires at most Lj + 1 ¥ -multiplications to 
compute the discrepancy and at most a further Lf + 1 < Lj for the updating when e« > 0. 
Ignoring subquadratic terms, this gives a worst-case upper bound of \\~\ ¥ -multiplications. 
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